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_ Abstract 

^ . We present a construction of a wavelet-type orthonormal basis for the space of radial 

' L^-functions in R'^ via the concept of a radial multiresolution analysis. The elements 

. of the basis are obtained from a single radial wavelet by usual dilations and generalized 

(N ■ 

translations. Hereby the generalized translation reveals the group convolution of radial 
functions in M'^. We provide a simple way to construct a radial scaling function and a 
' radial wavelet from an even classical scaling function on M. Furthermore, decomposition 

I ■ and reconstruction algorithms are formulated. 
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^ ■ 1 Introduction 

(N : 

I Standard approaches in multivariate wavelet analysis are based on the construction of mul- 

■ tiresolution analyses and wavelet bases from aflfine transformations of a finite set of basis 
, functions, called multi- wavelets (see e.g. |13| . [21). The translations are taken from a lattice 
I subgroup r of {W^, +) and the dilations are given by the integer powers of an expansive 
' matrix which leaves T invariant. The number of multi- wavelets needed to obtain a full basis 
. of L^(M'^) depends on the determinant of the dilation matrix and is in general larger than 

g ; 1 for d > 1. 

■ However, if one restricts to the analysis of radially symmetric functions only, it suggests 
K»" ■ itself to exploit this symmetry in the construction of corresponding wavelet transforms in 

^ I order to reduce the high amount of computational effort. A purely radial setting would 

^ ' for example naturally occur when separating variables in polar coordinates, and treating 

■ ■ the spherical and radial parts separately. There is a broad literature dealing with wavelet 

analysis and multiresolution on spheres, see e.g. PP, [HI and the references therein. In 
the radial case, it is not difficult to establish a continuous wavelet analysis based on the 
convolution structure of radial functions or measures instead of the usual translation in 
W^. Radial convolution structures are special cases (for half- integer indices) of so-called 
Bessel-Kingman hypergroups (see jl], J2]), and a continuous wavelet analysis can in fact be 
developed in this general setting, see e.g. ^H], and ^5]. Essentially the same concept 
is underlying the approach of Epperson and Frazier where radial wavelet expansions 
in are constructed which are based on sampling lattices with the spatial discretization 
determined by the positive zeros of related Bessel functions of the first kind. The spatial 
lattice is equidistant only in the special cases d = 1 and d = 3. This can be seen as an 
obstruction against a multiscale approach to radial wavelets in arbitrary dimensions. 
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As to the authors knowledge, radial multiresolution analyses have in fact not been con- 
sidered up to now, and there seems to be no general rigorous approach available for the 
construction of orthogonal radial wavelet bases in arbitrary dimension. This problem is 
closely related with the question to find a Poisson summation formula compatible with the 
Bessel-Kingman translation, which still remains open. We mention that the construction of 
j7j relies on the requirement that the involved radial wavelets are band-limited, i.e. their 
Fourier transforms have compact support. 

In the present paper, we construct radial multiscale analyses and orthogonal radial 
wavelet bases in M^. In dimension 3, the algebraic structure of the radial convolution 
allows to carry out the constructions along the same lines as in the well-known Euclidean 
setting. Hereby, the equidistance of the zeros of the corresponding Bessel function (which 
is simply a sinc-function) is of decisive importance. In order to motivate our approach 
and to make the intrinsic problems towards an extension to arbitrary dimensions more 
visible, we start the paper with a short account on the continuous wavelet transform for 
Bessel-Kingman hypergroups, and also present some material about Bessel frames. This 
in particular comprises the radial case in arbitrary dimension. The continuous transform 
is, up to normalization, the same as in ^H], while the Bessel frames are constructed in the 
spirit of the radial wavelet bases in jjj. In particular, in the spatial discretization the zeros 
of associated Bessel functions occur in a natural way, and the wavelets are band-limited. 

The radial analysis in is then also carried out in the setting of the corresponding 
Bessel-Kingman hypergroup H on [0, oo). Our concept of a radial multiresolution analysis 
(MRA) in is in fact that of a MRA for the L^-space L'^{H) of this hypergroup. The 
scale spaces {Vj)j^i C L^{H) are obtained by dyadic dilations from Vq, which in turn is 
spanned by equidistant hypergroup translates of a fixed "radial" scaling function (p G L^(H). 
It is characterized by a two-scale relation, but in contrast to classical MRAs, (j) itself is 
not contained in Vq, and the scale spaces are not shift-invariant (w.r.t. the hypergroup 
translation). Particular emphasis is put on the construction of orthogonal MRAs. Here 
periodicity arguments similar to the classical case are needed which would not be available 
in arbitrary dimensions. From a given orthogonal MRA we then derive an orthogonal 
wavelet basis for the underlying hypergroup. By construction, this "radial" basis has a 
direct interpretation as an orthogonal wavelet basis for the subspace of radial functions in 
L-^(R'^). We also provide a concise characterization of radial scaling functions in terms of 
even classical scaling functions on R. It in particular implies that in contrast to the classical 
case, there do not exist any real-valued orthonormal radial scaling functions with compact 
support. 

The paper is organized as follows: In Section 2, we recall basic facts from the analysis 
of radial functions in R*^, explain the corresponding radial hypergroup convolution struc- 
ture, and extend the setting to Bessel-Kingman hypergroups of arbitrary index. Section 3 
contains a short account on the continuous wavelet transform based on the Bessel-Kingman 
translation, as well as the construction of Bessel frames. In Section 4, radial multiresolution 
analyses in R^ and their scaling functions are introduced and discussed, while Section 5 is 
devoted to the construction of orthogonal radial wavelet bases. The connection between 
radial scaling functions in R^ and classical scaling functions on R is established in Section 
6. Finally, in Section 7 decomposition and reconstruction algorithms are discussed. 
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2 Radial analysis and Bessel- Kingman hypergroups 



Suppose F e L2(M'^) is radial, i.e. F{Ax) = F{x) a.e. for all A e SO{d). Then there is a 
unique / G L^{M.^,lO(1/2-i) such that F{x) = where | . | denotes the Euclidean norm 

on and for a > 0, the measure Ua on = [0, oo) is defined by 

dujair) = (2"r(a + l))~V2°^+^dr. 

Its normalization implies that ||-F||2 = H/lb, t^^/a-i' '^herc || . ||2 is taken with respect to the 
normalized Lebesgue measure {2Tx)~'^/'^dx on M"^. On L^(M_|_,a;a) the Hankel transform of 
index a is defined by 

/•oo 

r(A) = / U\r)f{r)dw^{r) 
Jo 

with the normalized Bessel function 

n=0 ^ ' 

There is a Plancherel Theorem for the Hankel transform, which states that / i— /" es- 
tablishes a self-inverse, isometric isomorphism of -L^(]R+, Wq,). If F G L^(R'^) is radial with 
F{x) = /(|a;|), then a short calculation shows that its Plancherel transform 

is again radial with J^{F){^) = /^'^/^"^Hl^D- This is due to the fact that 

Jd/2-i(kl)= / e-(^'f>da(0, (2.1) 

where da denotes the spherical surface measure normalized according to da{S'^~^) = 1. In 
contrast, the usual group translates x i-^ F{x + y) , y <E M.'^' , will no longer be radial (apart 
from trivial cases). However, we observe that the spherical means 



■F{x) := / F{x + rS,)da{S,), re 



of F are again radial. Moreover, ||MrF||2 < \\F\\2- Thus Mr induces a norm-decreasing 
linear mapping 

Tr : L\R+,ua/2-i) ^ L\R+,ua/2-i), Trfi\x\) := MrF{x), 

where / and F are related as above. Put a = d/2 — 1. Then a short calculation in polar 
coordinates gives 

Trfis) = Ca r /(Vr2 + s2-2rscos^) sin^'^ipdip with = ^/^"^w^ - (2-2) 
Jo 1 (a-l- 2ji (2) 

This defines a norm-decreasing generalized translation on L'^(R^,uOa) not only for a = 
d/2—1, but also for general a > —1/2. Having harmonic analysis in mind, we are thus lead 
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to introduce a corresponding measure algebra on M+r For r,s £ we define a probability 
measure 6r *a on M+ by 

Sr *a 6sU') ■■= Ca C f {^f + - 2rs COB ^) sin2°(^d(^, / G C,(M+). (2.3) 

(Cc(K_|-) denotes the space of continuous, compactly supported functions on M_|_). The con- 
volution (|2.3j) of point measures extends uniquely to a bilinear, commutative, associative and 
weakly continuous convolution on the space M{,(M_(_) of regular bounded Borel measures on 
M_|_. It is probability-preserving and makes -/Vfb(M_|_) a commutative Banach-*-algebra with 
respect to total variation norm, with neutral element 5q and the mapping fx >-^Jl as involu- 
tion. The pair (M+,*q,) is called the Bessel-Kingman hypergroup of index a. We write 
instead of M+ when putting emphasis on the specific convolution structure. Generally speak- 
ing, a hypergroup is a locally compact Hausdorff space together with a weakly continuous 
and probability preserving convolution of regular bounded Borel measures generalizing the 
measure algebra of a locally compact group; it also has a unit and an involution substituting 
the group inverse. In particular, every locally compact group is also a hypergroup. There 
is a well-established harmonic analysis for commutative hypergroups, which in the special 
case of reflects the harmonic analysis of radial functions (and measures) in M*^. We 

refer the reader to or |1] for a general background on hypergroups, including ^] for 
the Bessel-Kingman case. Let us mention only some aspects which are of importance in our 
context: The measure uJa is a Haar measure for H^, i.e. it satisfies 

/ Tsfduja= / fdu^ for ah / G C7,(M+). 
Jo Jo 

Up to a constant factor, uJa is the unique positive Radon measure on M_|_ with this property. 
Moreover, 

/ {Tsf)gduj^= / fiTsg)duJa for s G R+ (2.4) 
Jo Jo 

whenever both integrals exist. The Bessel functions satisfy the product formula 

5r *a Ss{ja) = ja{r)ja{s) for all r, s G M+; 

see |171 11.4]. For half- integers a, this is easily deduced from 1)2.1(1 . In fact, the functions 
f jai^i"), A G M+ are exactly those which are bounded and multiplicative with respect to 
*a- They constitute the so-called dual space of the hypergroup Ha- In this way, the Hankel 
transform / i— > /° on L^(M+,a;Q) can be interpreted as a Plancherel transform for Ha- For 
abbreviation, we put U'{Ha) := LP(M+,Wa) and we denote by (., .) and || . ||2 the scalar 
product and norm in L'^iHa)- It follows easily from ((TH) that for / G L'^{Ha), 

xTfiX) = ja{Xr)riX)- (2.5) 

3 Continuous wavelet transform and frames for Bessel-Kingman 
hypergroups 

In order to put the multiresolution approach in the following sections into a suitable frame- 
work, we continue with a short account on the continuous wavelet transform and wavelet 
frames for Bessel-Kingman hypergroups. This in particular includes a continuous radial 
wavelet transform and radial wavelet frames in arbitrary dimensions. 
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3.1 The continuous wavelet transform 

The following construction is essentially the same as in |16j , only with a different normaliza- 
tion of the dilation operators and the resulting wavelet transform. (In contrast to we 
choose dilations to be unitary, see below). We shall therefore be brief in our presentation, 
and refer the reader to [IBj Sect. 6. Ill] for further details. 

Let B{L'^{Ha)) denote the space of continuous linear operators on L'^{Ha)- Besides the 
translation operators Tr £ B{L?'{Ha)) introduced in Section 2, we consider the dilations 

Daf{r) «>0 
which are obviously unitary in L'^(Ha)- Notice also that 

(AJ)" = DyJ^. (3.1) 

We define 

TT-.H^x (0, oo) B{L\Ha)), 7r(r, a) := Tr-Da- 

It is easily checked that vr is continuous with respect to the weak operator topology on 
B{L'^{Ha)), c.f. Prop. 6.III.7]. 

Definition 3.1. A function g £ L^{Ha) is called admissible, if 
For abbreviation, put 

diJa{r,a) := ^2a+3 daduja{r). 
The following is a reformulation of |161 Thm 6.III.1] in terms of our notation: 
Theorem 3.2. (Plancherel Theorem) If g G L'^{Ha) is admissible, then for all f G L'^{Ha), 

[ K/,^(r,a)5)|'a!5,(r,a) = • ||/||2 . 

JHaX{0,oo) 

Polarization further implies for admissible gi, 52 and arbitrary /i, /2 S L^{Ha) the 
orthogonality relation 

dX 



A 



{fi,7r{r,a)gi){f2,7r{r,a)g2)duJa{r,a) = (/i,/2)- / 5f(A)?2°(-^) 
'Ha X (0,00) Jo 

Definition 3.3. Let g € L'^{Ha) be admissible. The mapping 

$g : L^{Ha) — > L^{Ha X (0,oo),5„), ^'g/(r,a) := {f,7r{r,a)g) 

is called the wavelet transform on Ha with analyzing wavelet g. For a = d/2 — \, it coincides 
with the continuous wavelet transform on of a radial function f with radial wavelet g, 

see rr^ . 

Some types of inversion formulas for this transform can be found in ^ 
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3.2 Bessel frames 



Let us now turn to possible discretizations. In order to obtain discrete versions of the usual 
wavelet transform on M, it is standard to use sampling lattices of the type {{nba^, a^), k,n ^ 
Z} with constants a > 1, 6 > 0. Here the discretization of the translation parameter is in 
accordance with the related group structure of M and is therefore (in general) not appropriate 
for radial wavelet analysis. Following Epperson and Frazier [Jj, we propose lattices where the 
discretization of the translation parameter involves the positive zeros Q < Va,i < '^a,2 < ■ ■ ■ 
of the Bessel function ja- By a result of McMahon, these are asymptotically distributed 
according to 

A standard lattice in Ha x (0, cxd) is given by 

{ {va,nha^ ,a^), A; G Z, n G N}, (a > 1, 6 > 0). 

The "almost orthogonal" radial wavelet expansions of Epperson and Frazier [TJ are based 
on this type of sampling lattice (with a = 2). In the following result, the discretization of 
the dilation parameter is still rather arbitrary. 

Theorem 3.4. LetQ be a countable subset of {0,oo) and g G L?'{Ha)- Assume that suppg°' 
is contained in [0, /] for some I > 0, i.e. = a.e. on (/, cxd), and that there exist constants 
A, 5 > such that 

-4 < ^ |?"(^A)p < B for almost all A G M+ . 

For n G Z and q € Q define "wavelets" gn,q G L?'{Ha) by 

gn,q ■■= ■ Tr^qDq{g) = 7r{rnq,q)g, 

where 



1 



-V, 



and M" 



2(i-")/2 J," 



t Vr(a + 1) I J„+i(z^„,„j| 

Then the set {gn,q '■ n. G N, g G Q} is a frame for L?'{Ha) with frame bounds A/^"^^ and 
BZ2"+2. This means that for f e L'^{Ha), 

Ai'-^'-\\m<EE\(9n,,j)\' <Bi'"^'-\\f\\i 

qeQneN 

Proof. The decisive point in the proof is the fact that the normalized Fourier-Bessel functions 

p^{X) := M„"j„(z/„,„A), nGN 

form an orthonormal basis of the Hilbert space Xa := ^^([0, 1], cjalfo,!]); see e.g. Erdelyi et 
al. Using (|TT|) and (^31) . we write 

?„",(A) = M^Ur^qX)Dyqg^iX) = p^{^-X)Dyqg"{X). 
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By the Plancherel theorem for the Hankel transform, we obtain 

where for the last identity, we used that the support of Di/i{'g°') is contained in [0,1]. 
Parseval's identity for now yields 

pi poo 

T.\(9n,gJ)\'= / \Dyi{g")D^/i{f-)fdcOa = / \g''{qX)\' \r{X)\' dua{X). 

_^,>,T JO Jo 



Hence, 



poo , , 



This implies the assertion. □ 

The second condition of this theorem is rather implicit. Following e.g. Bernier and 
Taylor [3, it is possible to obtain sufficient criteria which are easier to check by introducing 
the concept of "frame generators" . But in order to stay concise, we restrict ourselves to the 
most interesting special case of a standard lattice as defined above. Here Q = {a*^. A; G Z} 
with a > 1. 

Proposition 3.5. Suppose that g G Lp'^Ha) satisfies the following conditions: 

(i) has compact support which is contained in the open interval (0,oo). 

(ii) essinf\\g°'{\)\ : A G [a^jO""*"^]} > a for some n G Z and a > 0; 

(iii) r := ||5°||oo < oo. 

Then there exists a constant M > such that 

CT^ < J]] |5°(a''A)|2 < Mr^ for almost all A G M+ . (3.2) 

fcez 

Consequently, the set {gna'' '■ G N, A; G Z} is a frame for L^{Ha) with bounds cr^Z^"'''^ 
and Mr2/2"+2. 

Proof. We use the arguments of [3], sect. 4 in a simplified form which is adapted to our 
situation. Put F := [l,a). Then the intervals a^F = [a"',a^~^^), n G form a disjoint cover 
of (0, oo). As T := supp^" is compact in (0, oo), it is covered by finitely many of the a"F. 
This implies that 

M := sup ^{keZ: AG a^'^T} 

Ag(0,oo) 

is finite. By (iii), this gives the upper bound in (|3.'2j) . The lower bound follows from (ii) 
together with the fact that the a"'F cover (0, oo). □ 
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4 Radial Mult iresolut ion Analysis in 

Radial analysis in corresponds to the Bessel-Kingman hypergroup Ha with a = 1/2. For 
convenience we shall usually omit the subscript 1/2 and put 

duj{r) := duJi/2{r) = J -r'^dr, f := f^'^, j{r) := ji/2(r) = ^— . 

We further write H instead of -ffi/2 a.nd denote by ( . , . ) and || . ||2 the scalar product and 
norm in L?'{H) = L'^{H,ij), respectively. Notice that the Bessel function j is even on M. 
Hence it is natural to assume the Hankel transform 

/•oo 

fix) = / jiXr)f{r)du;ir) 
Jo 

of / G Lp'iH) to be continued to an even function on M as well. We shall always do this 
throughout the paper. We also mention that by a change of variables, the generalized 
translation (|2.2|) on H can be written in the simple form 

Trfis) = ^ f{t)tdt. (4.1) 

The non-negative zeros of the Bessel function j are given by 

tk := kiT, k e'N , 

and the normalized Fourier-Bessel functions 

Pk{r) := Mkjihr) with = 2^/\''/^k (4.2) 

form an orthonormal basis of the Hilbert space L^([0, l],u;|[o,i])- This is equivalent to the 
obvious fact that the functions 

s^{r) := (2/7r)^/^ rpfc(r) = \/2 sin(/c7rr), fcGN 

are an orthonormal basis for [0, 1] := -^^^([0, 1] , dr). It will be of importance in the following 
that the Sk are 2-periodic. 

Let us come to the definition of a radial multiresolution analysis (MRA) for W^, i.e. for 
the Bessel-Kingman hypergroup H. It is close to the well-known definition of Mallat Jl] 
for M. For convenience, we introduce the notation 

If / G L'^{H) then according to (^3]) . 

(M,rW/)-(A) = p,(A)/(A) = ilf-^fiX). (4.3) 

Definition 4.1. (Radial Multiresolution Analysis) A radial MRA /or M'^ is a sequence 
{Vj}j<^i of closed linear subspaces of L^{H) such that 

(1) Vj C Vj+i for all j G Z; 
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(2) nr=_oo = {0}; 

(3) Uj°l-oo is dense in L\H); 

(4) / G Vj if and only i//(2-) G Vj+i; 

(5) There exists a function 4> G such that 



:= {MkT^^U -.ken} 



is a Riesz basis ofVo, i.e. spanB^ is dense in Vq and there exist constants A,B>0 
such that 



< ||X]afcMfcTW(/<||i < B\\a\\l 



k=l 



for all a = {ak)km G ^^(N); here \\a\\2 = (Zl^i lafcP) 

The function (p in (5) is called a scaling function for the MRA {V^ }. We remark explicitly 
that in contrast to the classical case, (p itself is not contained in Vq and Vq is not shift 
invariant; in fact, if / G Vq then T^'^^f ^ Vq for all k. This will be shown in Corollary 14.61 
below. 

Our first aim is to determine an orthonormal basis for Vq from its Riesz basis, i.e. a 
function cj)* G L^{H) such that B^p* consititutes an orthonormal basis for Vq. For (p G Lp'{H) 
we define 

oo 

P^A) := 1^(^ + 2^)1'' 

n=— oo 

which is even and 2-periodic on M. 

Proposition 4.2. Let (j) G L'^{H) and A,B >Q. Then 

oo 

^llalli < II akMkT^^^ \\l < B\\a\\l for all a G /^(N) (4.4) 

k=l 

if and only if 

A<P^{\)< B for almost all A G M. (4.5) 
Proof. Let a G /^(N) be an arbitrary finite sequence. Define 

oo 

5:=^QfcSfc gL2[0,1]. (4.6) 

k=l 

We may regard a as an odd, 2-periodic function on M. By the Plancherel theorem for the 
Hankel transform and ()4.3() . 

oo oo oo 

^QfcMfcTW(/<||2 = ll^afcpfc^lli = / |^afcSfc(A)|2|^(A)|2dA 

7 1 / 1 7 1 



k=l k=l k=l 

|5(A)|2|(/>(A)|2dA = - I J5(A)|2p^(A)dA = / |5(A)|^P4A)dA. 



The Sfc forming an orthonormal basis of L^[0, 1], we have ||a!||2 = ||ct||£^2[o,i]- As the finite 
sequences form a dense subspace of /^(N), this implies the assertion. □ 
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With A = B = 1 we immediately obtain 

Corollary 4.3. For (j) G L^{H) the following statements are equivalent: 

(1) The set = {MkT^^^ : G N} is orthonormal in L'^{H); 

(2) P^ = l a.e. 

For 4> G L'^{H), put 

:= span_B<^ , 

the closure beine; taken in L^{H). The set B^ is a Riesz basis of if and only if there 
exist constants A,B > such that the equivalent conditions of Proposition are satisfied. 
This will be a standard requirement in the sequel, and we therefore introduce a separate 
notation: 

Definition 4.4. A function <p £ satisfies condition (RB) if B^ is a Riesz basis of 

As before, we shall often consider functions from L^[0, 1] as odd, 2-periodic functions on 
M. We therefore define 

S := {a : M — > C | a|[o_i] G ^^[0, 1], a{—x) = — a(x), a{x + 2) = a{x) for almost all x}. 

is a Hilbert space with norm || . ||l2[o_i]. 

Lemma 4.5. Let cp G satisfy (RB). Then for f G Lp'{H) we have the equivalence 

feV^^ /(A) = ^0(A) with (3eS. 

The function / G corresponding to (5 = Yl'kLi '^kSk £ S with {ak)keN G /^(N) is given by 

Proof. By (|4.3() we have 

^^(A) = {^f\MkT(''UrW. 

The translates M^T^^^cj) form a Riesz basis of V^, and hence the (Mf^T^^^cf))^ are a Riesz 
basis of V(j,. As {sk)km is an orthonormal basis of S, this implies the assertion. □ 

Lemma 14.51 is of particular interest when (p is the scaling function of a MRA {Vj}. Then 
Vq = V^, and we easily deduce the previously mentioned lack of shift-invariance: 

Corollary 4.6. Let {Vj}j(zz be a radial MRA. Then f e Vq implies that T'^'^^f = TtJ i Vq 
for all G N. Similarly, f & Vj implies that T2-jt^ f ^ Vj for all A; G N. 

Proof. After rescaling it is enoug h to consider Vq. Recall that Mk{T^''^f )^ = pkf. But if 
f3 £ S, then pkP ^ S for all k, because periodicity is lost. The characterization of = Vq 
according to the previous Lemma thus shows that for / G Vq, T^^'^f ^ Vq. □ 

In the situation of the Lemma, we can easily determine an orthonormal basis of by 
renormalization: 
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Theorem 4.7. (Orthogonalization) Suppose (j) G L?'{H) satisfies condition (RB). Define 
(jf G LP'{H) by its Hankel transform 




Then B^* = {M^T^^^^cp* : G N} forms an orthonormal basis ofV^ = V^* . 

If (/> is a scaling function of a MRA {V^ }, then = V^* = Vq, and we call (j)* an orthonormal 
scaling function for {Vj}. 

Proof. By definition of 0* we have P,^* = 1 a.e. and hence i?,^* is orthonormal according 
to Corollary 14.31 It remains to prove that V^* = V^. For this, we have to verify that 
MfcT^'^V* G V4, and MkT^^^ G V^. for all /c G N. Employing Lemma ESI relation 
and finally the relation Sfc(r) = (7r/2)-'^/^ rpfc(r), one obtains that the above conditions are 
equivalent to 

G S, Sky^ G S for all A; G N. 

But these conditions are obviously satisfied by our assumption on . □ 

Let us return to our definition of a radial MRA for M'^. Suppose we start with a function 
(j) G L^{H) satisfying condition (RB) with Riesz constants A,B>0. Define corresponding 
scale spaces {Vj}ji^z by 

where the dilation operator Da G B{L^{H)) is defined as in Section 3. Then in particular, 
the Vj satisfy axiom (4) of Definition 14.11 Put further 

0,-fc(r) := D2-, {MkT^''U){r) = Mk{T^^^ (t>){2^ r) j£Z,kGn. (4.8) 

Then {(pj^k, = {4>o,k-, 4'o,i) for all j, k, I. Thus the {(/'j,^, A: G N} form a Riesz basis of Vj, 
with the same Riesz constants A, B as for j = 0. In particular, 

Vj = span{4>j^k,k G N}. 

Moreover, if B^ = {(po^k '■ A; G N} is an orthonormal basis for Vq, then {(pj^k, A; G N} is an 
orthonormal basis of Vj. 

Recall now axiom (1) of Def. 14.11 which requires that the Vj are nested. Similar as in 
the classical case, this condition can be reformulated in terms of a two-scale relation for (p: 

Proposition 4.8. For cp with (RB) and {Vj}j,^z as above, the following statements are 
equivalent: 

(1) Vj C Vj+i for all j G Z 

(2) F-i C 

(3) There exists a function 7 G such that 

sin(27rA)^(2A) = 7(A)^(A). (4.9) 
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In this case, the coefficients {hk)keN £ ^^(N) in the two-scale relation 

oo 

<^-i,i = ^hk(^o,k (4.10) 

k=l 

are the coefficients in the Fourier sine series 0/7 G S: 

1 °° 



00 

i 

7 



2 

k=l 

Proof. Rescaling by the factor 2^ shows that (1) and (2) are equivalent. For (2), we need at 
least (/>-i,i G Vq. According to Lemma 14.51 this is equivalent to the existence of a function 
(3 ^ S such that for almost all A, 

\/8Mij(2vrA)^(2A) = ^i(A) = ^0(A). (4.11) 

Moreover, if has the expansion (|4.1flj) . then (3 = {tt /2Y^^ Yl^=i ^kSk- In turn, equation 
()4.11() is equivalent to relation 1)4. 9() with 7 = (87r)-i/4/?. This gives the stated connection 
between 7 and It remains to show that € Vq (or equivalently, relation (|4.9|) ) 

already implies that 0_i^fc G Vq for all A; G N. Similar as above, the latter is equivalent to 

sin(2A;7rA)^(2A) = A;7fc(A)^(A) (4.12) 

with 7^, G S. The relation between 7^ and 0-i fc is now given by 



7fc = -^^^f^^k, 4>-i,k = ^h\^''4>o,i- 



1=1 1=1 
Comparison of (|4.9|) with (|4.12l) yields 

7fe(A) = iW^^^^^^, = 7(A)^fc-i(cos2vrA) (4.13) 
k sm(27rAj 

where 

sin(A; + l)t 
Uk{x) = , . . , x = cost 
(k + 1) sm t 

denotes the fe-th Chebychev polynomial of the second kind, normalized such that Uk{l) = 1. 
Thus given 7 G 5", we define 

7fc(A) := 7(A)[/fc_i(cos27rA). 
As t/fc-i is bounded on [—1, 1], 7^ is contained in S as well, and hence (p-i^k ^ ^o- D 
Let us now consider the remaining axioms (2) and (3) of a radial multiresolution analysis. 
Theorem 4.9. Let (p G Lp'{H) satisfy condition (RB) and assume that the scale spaces 



Vj = span{(pj^k : fc G N}, i G Z 

satisfy V-i C Vq. Suppose further that \(f>\ is continuous in 0. Then {Vj}j^z is a radial 
MRA if and only if 0(0) 7^ 0. Moreover, <j) is an orthonormal scaling function if and only 

if \m\ = 1- 



12 



We remark that continuity of </) in (even on is for example guaranteed if (/> G 
L'^iH)nL^H). 

Proof. We have to check axioms (2) and (3). This may be done by shght modifications of 
standard arguments in the affine case. The condition on (p in will be needed only for (3). 
We define an orthonormal scaling function (f)* according to Theorem 14.71 The orthogonal 
projection Pj of L'^{H) onto Vj is then given by 



klVj,k 



where the (jf- ^ are defined as in (|4.8j) . For (2), we need to show that limj^_oo H-Pj/lb = 
for all / G LP'{H). Since functions with compact support are dense in L'^{H), we may assume 
that supp/ is contained in a compact interval [0, ii]. Parseval's equation then implies 



k=l k=l 

fR °° (■2m °° 
<\\f\\l- J2\^lk{r)\'du;{r) = \\fg- ^\cPl,{r)\'dco{r). 
■^0 k=l fc=l 



Using the explicit formula (|4.1j) for the hypergroup translation in H we further deduce 

-(1/4 fkn+r 2 

_ _ (l)*{t)tdt 

k=l ' k=l fc=l' J\kn-r\ 



00 00 00 /r) xl/4 fkn+r 



Now assume that j is sufficiently small so that 2^R < tt/2. Then for r E [0,2-^i?], the 
integration domains [kn — r, /cvr + r] do not overlap and we obtain 



El'^S..(OP < 7E/ W{t)\H'dt<'' 



* l|2 
— IIV^ 2 
r 



k=l k=l' 

with suitable constants C, C" > independent of j. Hence for j sufficiently small, 



\Pjf\\l < C" / -dcoir), 
Jo 



which tends to as j ^ — oo. This proves (2). 

As to (3), suppose first that 0(0) / and let /i G ( U^_^ Vj)^ , i.e. Pjh = for all j £ Z. 
We claim that h = 0. Indeed, for e > there exists a function / G L?'{H) such that the 
support of its Hankel transform / is compact and ||/ — /i||2 ^ £• This implies 

Il^':,-/ll2 = \\PjU-h)h <e for ah j G Z. 
By the Riesz basis assumption on (j), we further have 

oo oo 

AY,\{fA,,k)? < WPjfWl < BY,\{f,^j,k)\\ (4.14) 

k=l k=l 
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c.f. Lemma 2.7 in jl8j . Further, if supp/ G [0, i?] then 



R 



{f,4>j,k) = {f,$,,k)= / /(A)p?(A)0(2-U)d^(A) 

Jo 

where 

Note that the functions {pl^\ A; € N} form an orthonormal basis of i^([0, 2-^], u;|[o,2j]) =• -^j- 
Suppose now that j is sufficiently large, i.e. 2^ > R. Then 



{f,<t>j,k) = (/?(2-^-),p1''^)x,. 
Thus by Parseval's equation for Xj, 

oo f-R 

EK/'-^^^)!' = l|/?(2"'-)|li. = / |/(A)n0(2-^A)pda;(A). 
k=i ' 

As is assumed to be continuous in 0, the functions A \4>{2^^X)\ converge to the 
constant \4>iG)\ > uniformly on [0,R] as j oo. Hence 

e > lim ||P,-/||2 > VA\m\\\fh > ^\m\i\\h\\2-e). 

As e is arbitrarily small, this shows that /i = and hence axiom (3) is satisfied. Vice versa, 
axiom (3) implies that 

lim Pjf = f for all / G L^{H). 

If / is compactly supported, then the same calculation as above shows that 

lim ||P,/||2 < VB\m\\\f\\2 

j->oo 

which enforces ^(0) 7^ 0. Finally, the <j)j^k , k £'N are orthonormal if and only if it is possible 
to choose A = B = 1 in (|4.14|) . In this case we obtain (for / as just before) 

lim ||P,-/||2 = |0(O)|||/||2 = 1^(0)1 II/II2. 

Thus (3) is satisfied exactly if 10(0)1 = 1. □ 
Let us now write the two scale relation 1)4. 9() in a slightly different form, namely 

0(2A) = G(A)0(A) (4.15) 

with 

G(A):^ 



sin(27rA) 

The filter function G is obviously 2-periodic and even. Whenever its restriction to [0, 1] is 
contained in -L^[0, 1] it can be expanded as a cosine series, 

00 

G(A) = V^y~^ ffn cos(n7rA). 

n=0 

As for a classical multiresolution analysis one proves the following. 
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Lemma 4.10. Suppose that (j) G L'^{H) is an orthonormal scaling function of a radial MRA. 
Then the associated filter function G satisfies 

|G(A)|2 + |G(A + 1)|2 = 1 a.e.. (4.16) 

// in addition (j) £ (H) then ^4-16^ holds pointwise and 

G(0) = 1, G{1) = 

which implies 

oo oo 

V2^5n = 1, V2^(-ir5n = 0. 
n=0 n=0 

Proof. In view of Corollary 14.31 we have 

oo oo 

1= Yl l^(A + 2n)|2= |G(A/2 + n)|2|0(A/2 + n)|2 

n=— OO n=—co 

oo oo 

= |G(A/2)|2^ |^(A/2 + 2n)|2 + |G(A/2 + l)|2 ^ |0(A/2 + 2n + 1)|2 

n=— oo ?i=— oo 

= |G(A/2)p + |G(A/2 + l)p 

almost everywhere. If G (H) then (p is continuous and 0(0) 7^ by Theorem 14.91 
Hence G(0) = 1 by (|4.15j) and G(l) = is an immediate consequence of (j4.16|) . □ 

5 Orthogonal radial wavelets 

In this section we construct wavelets for a given radial MRA {Vj}ji^z in with orthonor- 
mal scaling function and filter function G. As usual, the wavelet space Wj is defined as 
the orthogonal complement of Vj in Vj+i, 

Thus LP'^H) decomposes as an orthogonal Hilbert sum 



00 

L\H) = W,. 



j=-oo 

Recall the definition of 5 in Section 4 and the characterization of Vq = K/> according to 
Lemma 14.51 

feVo^ /(A) = ^ 0(A) with pes. (5.1) 

Define 

Sq := {a £ S : a(A + 1) = — a(A) for almost all A} 

which is a closed subspace of S w.r.t. || -111,2 [g^ij. Then W^i = Vq V-i is characterized as 
follows: 
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Proposition 5.1. (i) Let f E L^{H). Then 



f E W_i ^ /(A) = G{X + 1) ^(A) for some a G So- 

A 

(ii) The mapping Sq W-i ^ fa with 

UX) := (2vr)i/4^G(A + l)0(A) 

is an isometric isomorphism. 

Proof. For (i), notice first that the Hankel transform is a unitary isomorphism of L'^{H), so 
W-i = VqQV-i. Rescaling of by the factor 2 and relation (I^TT^ imply that h £ L'^{H) 
is contained in V-i if and only if there exists some f3 G S such that 

MA)=M?,2A) = Mg„)J(„. 

Thus P G S corresponds to / S W^i according to ()5.1|) if and only if 



/?(2A) ^^^^ ^|^(^)|2^^(^) ^ for all /3 G 5. 



A ' ' A 

Up to a constant factor, the integral on the left equals 



/?(2A)G(A) /5(A) |(/)(A)p dA = V / /3(2A) /3(A + 2n)G(A + 2n) |0(A + 2n)|2dA 
+ V /" /3(2A) /3(A + 2n + 1)G(A + 2n + 1) |^(A + 2n + l)|2dA 



n&L' 
-1 



/3(2A)(/5(A) G(A) + /3(A + 1) G(A + \))dX, 

where we used the periodicity and symmetry properties of I3.,I3^G as well as Corollary 14.31 
Since /? G S" is arbitrary, we conclude that the vectors (/3(A), /3(A+1))"^ and (G(A), G(A+1))"^ 
must be orthogonal in for almost all A. This means that 



for some function a : [0, 1] C. Thanks to the boundedness of G, a belongs to -L^[0, 1]. 
Since [3 and G are 2-periodic on R, (3 is odd and G is even, an extension of a to M must be 
2-periodic and odd; hence a G 5 with /3(A) = a(A)G(A + 1). Using the 2-periodicity of (3 
and G in the second component of 1)5. 2|1 . we further deduce that /3(A) = — a(A + 1)G(A + 1) 
and therefore a(A + l) = — a(A) apart from the zero-set of G(A-|-1). Thus / is of the claimed 
form. Conversely, if a G 5*0, then /3(A) := a(A)G(A + 1) G S*, and (|5.2|) is satisfied. 
For the proof of (ii), we calculate 

Wf^Wl = wig =V2^jJ|"M^|G(A + l)|2|^(A)pd^(A) 

|a(A)|2|G(A + l)|2|^(A)|2dA = j' |a(A)|2|G(A + l)|2dA, 



16 



where we used that a and G are 2-periodic and </> is an orthonormal scahng function. By 
assumption on q, we have a(A — 1) = a(A + 1) = — a(A). Thus by Lemma I4.1U1 

1 ,\ 1-1 

|a(A)p|G(A + l)|2dA = / |a(A-l)|2|G(A)|2dA + / \a{X)\'^\G{X + l)\'^dX 



1 



1 

|a(A)pdA = \\a\\l. 



This proves (ii). □ 

It is now easy to obtain an orthonormal basis of W-i. Recall that the Sfc, A; G N form an 
orthonormal basis of S. Moreover, the S2k~i, k £ N are an orthonormal basis of 5o. Thus 
by the previous result, the functions 

constitute an orthonormal basis of W-i. Define ip G L^{H) by 



V'(2A) = G(A + 1)(/>(A). (5.3) 

Then in view of (|4.3() . 

To obtain an orthonormal basis of Wq, we just have to rescale. Extending the notation 
r('^) := to r G N/2 and using the relation DaT^ = TaxDa, we obtain that an orthonormal 
basis of Wq is given by the functions 

V'. := Dy,.h = M^T^i^^-^m^, k G N. 

We call ij) a (basic) wavelet for the radial mutiresolution iVj). 
Definition 5.2. For j G Z and A; G N, define the "radial" wavelets 

V.,. ,(.) := D,-Mr) = 8^V2^!^|zir((2'=-i)/2)^(2i,). 

We have proven: 

Theorem 5.3. (i) For each j G Z, the set {ipj^k : A; G N} constitutes an orthonormal 
basis of Wj . 

(ii) The set {ipj,k : j G /c G N} is an orthonormal wavelet basis of L'^{H). 
Corollary 5.4. The functions 

^j^k{x) := ^Pj,k{\x\), xeR^j ez,ken 

form an orthonormal basis for the closed subspace L^^^(M^) := {/ G L^(M'^) : / radial} of 
radial functions in L?' 
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6 Construction of radial scaling functions and wavelets 



Yet, we do not have handsome criteria in order to decide whether a given function (j) G Lp'{H) 
is a radial scaUng function, i.e. a scaUng function for a radial MRA. The analogy of our 
constructions to those on the group (M, +), however, leads to the following close relationship: 

Theorem 6.1. Suppose (p^ is a classical scaling function on M which is even and such that 
its (classical) Fourier transfrom is continuous in and satisfies J-{(j)u) G L^{H). 

Define 4> € -^^(^) ''^^'^ ^^-s Hankel transform, 



</.(A) := V2^.^(0k)(^A). (6.1) 

Then (j) is a radial scaling function. 

Conversely, if (j) is a scaling function for a radial MRA such that <j) is continuous in 
0, then (j) £ L^(M) and the function cp^ defined by id.l]) (where (p is extended to an even 
function on M) is a classical scaling function on M. 

Moreover, (p is an orthonormal radial scaling function if and only if (p^ is an orthonormal 
classical scaling function. 

Proof. Let us start with the first assertion. As (p^ is a classical scaling function, we have by 
eq. (5.3.2) in 

fcez 

with suitable constants {) < A < B < oo. Moreover, c^r is orthonormal if and only if 
A = B = 1. Since (p^ is assumed to be even, definition ()6.1() is compatible with the even 
extension of <p. By Proposition 14.21 the set {Mj^T^^^(p, k E N} forms a Riesz basis for Vq = 
span{i?(^} which is an orthonormal basis if and only if (p^ is orthonormal; c.f. Corollary 14.31 
Moreover, by eq. (5.3.18) of there exists a 27r-periodic function rriQ £ L^([— 7r,7r]) such 
that J^{<pM.){(,) = "T'o(C/2).^((/'r)(^/2), and mo is necessarily even in our case. Hence, with 
7(A) := mo(A7r) sin(27rA) which clearly is contained in S, we have sin(27rA)(/)(2A) = 'y{X)(p{X). 
This is exactly the radial two-scale equation ()4.9|) . As (p is continuous in 0, the condition 
0(0) 7^ of Theorem 14.91 is automatically satisfied (see e.g. Remark 3 on p. 144 in 6 ), and 
thus we finally obtain that (p is a radial scaling function. 

For the converse part notice first that continuity of i?!> in already implies that p £ (M) . 
We further proceed similar as before, using Proposition 5.3.1 and Proposition 5.3.2 in |S] 
and the corresponding results of the present paper. □ 

This theorem supplies a variety of radial scaling functions since there are many classical 
scaling functions on R which satisfy the assumptions of the theorem. However, as to or- 
thonormal radial scaling functions with compact support, a famous theorem of Daubechies 
implies the following negative result. 

Corollary 6.2. There do not exist any real-valued orthonormal radial scaling functions with 
compact support. 

Proof. The proof of Theorem 8.1.4 in shows that an even, real-valued and compactly 
supported scaling function is necessarily the Haar function X[-i/2,i/2]) the characteristic 
function of the interval [—1/2,1/2]. However, its Fourier transform J^{x[-i/2,i/2]){0 = 

./ ig not contained in L^{H). □ 
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If (f) correponds to an even classical scaling function (f)^ according to Theorem 16 .11 then 
the hypergroup translates 0o,A: = ^kTk-w'P may be expressed accordingt to the formula 

Mx) = ^^-i^ (<^k(^ -k)- <Am(^ + k)) . (6.2) 
In fact, by the Plancherel theorem for the Hankel transform and eq. (|2.5|1 . 

) -1 poo 

j{\r)j{\s)4>{\)du{\) = -j= / ^(A) sin(sA) sin(rA)dA 

V27r • rs J-oo 



for all r, s G M+. Here (p is as usual extended to an even function on M. Using relation H6.1|) 
as well as basic trigonometric identities and the Plancherel theorem for the classical Fourier 
transform, we can write 

Tr4){s) = — j= / (^(A)(cos A(r — s) — COS A(r + s))(iA 



2\/27f • rs 



1 /, (r — s\ , (r + s 



27r • rs V V vr / V tt 

This implies (|6.2() . 

As an example, we consider the radial analogue of the Shannon wavelets. We define the 
scaling function via its Hankel transform, 

2i\\ (W \ /2" sin(a;) -xcos(x) 
0(A) = X[o,i](A), 4>{x) = ^3 . 

Constructing the associated basic wavelet according to formula (|5.l-i|) yields (after a short 
calculation) 

; /2 sin(2x) — sin(x) — 2x cos(22;) + X cos(x) 
V'(A) = X[i,2](A), ^{x) = ^3 . 

The translates of the scaling function and the wavelet turn out to be 
, , 1 /sin(x - /cvr) sin(x + /cvr) \ , 

, , ^ 1 ( ( 2A:-1 \ / 2A:-1 \\ 
V'0,fc(2;) = ,^ ,1/, p X — TT -p\x^ TT , A; G N, 



with ^{x) 



sin(2x) — sin(a;) 



X 



7 Algorithms 



For the use of our radial multiresolution in applications we need to formulate decomposition 
and reconstruction algorithms. The first step in such an algorithm consists of projecting 
the function / into a scale space Vj for some suitable j. We obtain a representation 

oo 
k=l 
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So from now on we assume that we have given a function / E V^- in terms of its coefficients 

(i) 

cj^ . The decomposition algorithm consists of decomposing / into Vj-i and Wj-i, i.e. of 

(j—i) (7 — 1) 

calculating the coefficients and dj^ in the representation 



00 



k=l k=l 

(Such a representation exists, since by construction {^j-i^fc, '>Pj-i,k '■ A; G N} is also a basis 
of Vj.) A reconstruction algorithm determines the coefficients cj^^ when / is given in terms 

of cjf "^^ and djf "^^ keN. 

We still assume that (j) is an orthonormal scaling function (and ip is hence an orthonormal 
wavelet). Let 

■= {4>i,k,4>0,l) = {4'j,k,4'j~i,l), := (01,^,-00,/) = {4'j,k,il^j-i,l)- 

By using Hilbert space techniques - in particular Parseval's equation - we obtain analogously 
as in standard wavelet theory the decomposition formulae 

00 00 

^£ ~ Z^^fc ' ~ 2-^^k 'i ' 

fc=i fc=i 

and the reconstruction formula 

00 00 

£=1 i=\ 

It turns out that the coefficients qf"^ and r^*"^ are determined in terms of the numbers gn in 
the cosine expansion of G, i.e. the coefficients in 

00 

G(X) = \/2 gn cos(n7r A) . 

n=0 

Theorem 7.1. For £,k eN it holds 



(k) 
1i 



gk-2e - 9 2£+k for 21 < k, 
250 - 9M for 2£ = k, 



, g2e-k - g2i+k for 2£ > k, 

{-lf~^{gk^2i+i - 9k+2i-i) for 21 -Kk, 
2go - 9U-2 for 2i-l = k, 

{-l)''^\g2e^i-k - 92£-i+k) for 2^ - 1 > A;. 
Proof. Using the Plancherel theorem, relation 1)4. 15(1 and Corollary 1)4. 3|) we obtain 

={4'i,k,4'o,e) = {(i>i,k,(i>o/) 

"1/2 / n,(\l'}\n„(\l9'\n(\l9\\^(\l9'\\'^, 



Pk{\/2)p,{\/2)G{\/2)\(t){\/2)\'du:{\) 
■V2 Sk{\)si{2\)G{X) \^{\ + 2n)\^dX 

•^^ n=-oo 
00 „i 

-Ay^g^ / sin(A;A7r) sin(2£A7r) cos(n7rA)(iA. (7.1) 

n=0 
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An easy calculation using trigonometric identities shows 



sin(/i:A7r) sin(tA7r) cos(n7rA)dA 



t-k\ 

3(25o,t-fc 



^n,t+k) for n > 0, 
(5o,t+fc) for n = 0. 

(fe) 



(7.2) 



Setting t = 21 and inserting into (|7.1() yields the assertion for . We proceed similarly for 



rf =(</'i,fc,V'o/) = ^ / Sfc(A)s2£-i(A)G(A + l)dA 

JO 

sin(A;7rA) sin((2£ - l)7rA) cos(n7r(A + l))dA 

oo „i 

4V5„(-1)" / sin(/fc7rA)sin((2£- l)7rA)cos(n7rA)(iA. 



oo 

n=0 

oo 



Setting t = 2£ — 1 in ()7.2() and inserting into the last expression gives the result for r'l^ . □ 

Let us consider the case where only finitely many coefficients are different from zero. 
Although this is not possible for real-valued orthonormal scaling functions this assumption 
makes it easier to compare the radial wavelet algorithm with the classical one. Of course, 
in applications one can only handle finitely many coefficients anyway. So let us assume 
suppg C [0,A^], i.e. = for k ^ {0, ...,A^}. Elementary considerations show the 
following. Leaving k fixed yields 



(k) 
U 

I 



for£^[^,^] 

for £ ^ [1, ^] 

for £ ^ [^^, 

for ^ ^ [1, ^l±k±l] 



ifk>N, 
iik<N, 
i{k> N + 1, 
iik<N + l. 



If £ is fixed then 



(fc) 
U 

(k) 

u 

Sk) 

I 



With 



for /fc ^ [2£ -N,2i + N] 

iov k ^[l,N + 21] 

for /c ^ [2^ - 1 - iV, 2i - 

for A; ^ [1,2^ - 1 + iV] 

g\^ for 1 < \k\ 
2gQ for k = 0, 
otherwise. 



if 2i > N, 
if 2i < N, 
l + N] if 2^ - 1 > iV, 
if 2£ - 1 < TV. 

<N-1, 



it holds G{X) = ■^'Ylik=-N^k(^^^^ ■ Because of the conditions on it, G is also the filter 
function for an ordinary multiresolution analysis on M with coefficients /i^. Now, if 2^ > 
resp. 2^ — 1 > A^ then it is easy to see that 



=hk for A; = -A^, . . . , N, 
^fl-i+k) for fc = -A^,. 



21 



Similarly, ii k > N then 

) = h^e, = {-lf-+'hk+i-2e. 

Hence, for 2i > N + 1 the decomposition formulae become 

oo oo 
k=l k=l 

and for A; > iV + 1 the reconstruction formula is 

oo oo 
£=1 £=1 

These formulae are well-known. Indeed, they are the decomposition and reconstruction 
formulae of the classical discrete wavelet transform. So our approach leads to the classical 
algorithm if we are far enough away from the origin. If we are close to the origin we have 
derived an algorithm to handle the boundary point 0. 
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